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Numerical experiments on dynamo action in sheared and rotating 
turbulence 
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Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical 
turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field 
at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate 
S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed 
shear, 7 oc 5*, while the characteristic spatial scale of the field is Ib oc S~'^^^. The effect is quite general: earlier results for 
the nonrotating case by Yousef et al. (2008) are extended to shearing boxes with Keplerian rotation; it is also shown that the 
shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic 
nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of 
magnetic fields in astrophysical systems. 
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1 Introduction 

Turbulence is generally considered to play a fundamen- 
tal role in the generation and maintenance of mag- 
netic fields found in a wide range of astrophysical 
systems. It is a numerically well established prop- 
erty of turbulence to amplify magnetic fluctuations at 
the same or smaller scales than the scales of the 



turbulent mo t ions (IHaugen. Brandenburg & Doblei 



Iskakov et al. 2007 



Meneguzzi. Frisch & Pouqueu 



2004 



1981 



Schekochihinetal.ll2004l2007h . This type of dynamo is 
known as small-scale, or fluctuation, dynamo. It is belieyed 
to be a uniyersal property of turbulent systems and, at least 
in the case of large magnetic Prandtl numbers, a simple the- 
oretical picture exists of the field amplification yia random 
stretching by turbulent motions (IMoffatt & Saffmanlll964l: 
Zel'doyich et al. 19841: see also ISchekochihin & Cowley 



2007: Schekochihin et al.1 120041) . It often turns out to be 



more difficult to either establish or explain the generation 
of magnetic fields at much larger scales than the turbulence 
scale. Such fields are obseryed or belieyed to exist in many 
astrophysical bodies (stars, galaxies, disks), so the question 
of their origin is an important theoretical challenge. 



Mechanisms for the generation of such large-scale 
fields are known as large-scale, or mean-field, dynamos. 
A motley of such mean-field dy namos has been studied 
in the literature (e.g., iBra ndenbur g & SubramanianI 12005 



^ Corresponding author: tarek@pvv.org 



iK rause & Ra"dyil980l: lMoffattiil978[]Radler & Rheinhardtl 



2007h . The precise way in which they operate often seems 
to be system dependent and introducing eyer more real- 
istic features into one's theoretical model produces eyer 
more complicated behayiour. While such modelling is nec- 
essary for quantitatiye understanding, it is quite interest- 
ing to ask what are the generic ingredients required to pro- 
duce large-scale fields. One such ingredient appears to be 
the presence of net kinetic helicity in the system: in many 
mean-field dynamos, the key generation mechanism is 
the so-called a effect dSteenbeck. Krause & Radler 1966), 
whereby an assembly of non-mirror- symmetric yelocity 
fluctuations haying a nonzero net helicity are responsible 
for magnetic-field generation. The existence of a mean-field 
dynamo in helically forced turbulence is well established 
numerically CBrandenbur glbOO ll: iBrandenburg et al.ll2008bl: 
Maron & Bla ckman 2002). Howeyer, the requirement of net 
helicity may somewhat limit the applicability of the a ef- 
fect. It is also far from certain that the direct link between 
kinetic helicity and mean-field generation yia the a effect 
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found in the model case where the helicity is injected by 
the random forcing, carries over to the cases where the he- 
lici ty arises naturally (e.g., in a rotating convective layer; 



Cattaneo & Hughesll2006l) . Therefore, there is a strong 



see 

motivation for seeking alternative mean-field dynamo pro- 
cesses driven by nonhelical turbulence. 

It is clear that a nonhelical turbulence by itself can- 
not make large-scale fields. In recent years, many au- 
thors have argued that large-scale magnetic fields can be 
generated by nonhelical velocity fluctuations when acted 
upon by a large-scale shear: theoretical paradigms put 
forward to support such a dynamo, which we refer to 

as th e shear dynamo, have included the sh ear-current ef- 
I • ' — -7; 1 1 — • — I I 1 

feet (iRogachev skii & Kleeorin 2003, 200 4D, the st ochas- 
tic a effec t fBran denburg et al. 2008a; Ipedotoy 
Fedotov. B ashkirtse va & Rvashkol 120061 IProctoJ 



nn 



2003 



T 



120071: 



Silant'e v 2000; Vishniac & Brandenburg 1997; see, how- 
ever, iKleeorin & Rogachevskii 2008), negative-di ff usivitv 



type theories dRiidiger & Urpinll200ll: IUrpinlll999allb[ 12002 , 
20061), shear amplification of the fluctua tion-dynamo- 
generated small-scale fields (lBlackmanlll998b . While there 
is not yet agreement between theoreticians about the valid- 
ity or areas of practical applicability of these models, it is 
clear that they are addressing a fundamental issue. Indeed, 
shear is an extremely generic property of astrophysical sys- 
tems, so the idea of a shear dynamo gives us a particularly 
attractive scenario for ubiquitous generation of large-scale 
magnetic fields. 

Until recently, a numerical demonstration of this type of 
dynamo re mained elusive. O riginally motivated by the pre- 
dictions of Ro gachevskii & Kleeorin (2003, 2004), we have 
previously performed numerical simulations of nonhelical 
turbulence with a superimposed linear shear and demon- 



strate d the existence of the shear dynamo (lYousef et al. 



2008|) . Theoretical understanding of these numerical results 
is still poor. More work and, we believe, more information 
gathered from numerical experiments are needed in order 
to make progress towards understanding the properties of 
this dynamo and the underlying physical processes that pro- 
duce it. This paper , which is an extension of the work by 
Yousef et al.l (l2008k . aims at presenting a collection of new 
numerical results regarding the existence and behaviour of 
the shear dynamo in various regimes. We focus on three dif- 
ferent cases of astrophysical interest, namely shear dynamo 
in the presence of forced nonhelical turbulence and a linear 
velocity shear, shear dynamo in the presence of forced non- 
helical turbulence and Keplerian differential rotation, and 
finally shear dynamo in the presence of forced nonhelical 
turbulence, Keplerian differential rotation and a small-scale 
fluctuation dynamo. In order to study the effects of shear 
and rotation, we adopt a local rotating shearing sheet model. 
This model and the corresponding numerical set-up are pre- 
sented in Sect. [2l In Sect. [3l we consider the case of linear 
shear without rotation. Results for the Keplerian regime are 
presented in Sect. HI Section [5] describes some preliminary 
results on the shear dynamo in the presence of small scale 



magnetic fluctuations generated by the fluctuation dynamo. 
A short discussion concludes the paper (Sect.©. 

2 Model and numerical set up 

We consider differentially rotating flows that can locally be 
described in terms of a background shear flow U = Sxy 
rotating uniformly with a rotation rate Viz. In the frame- 
work of incompressible magnetohydrodynamics (MHD), 
solenoidal velocity and magnetic field perturbations to this 
flow evolve according to the following equations written in 
the rotating frame: 

Vr) 

(1) 



du 
dt 


= —u 


j,{S ^2n)y ^2nuyx - 


Vp 
p 




^B 


VB^uV^u^f, 




dB 
dt 


= B, 


,Sy^B-Vu^r]\/^B, 





(2) 

where d/dt = d/dt + (t/ + n) • V, n is the velocity de- 
viation from the background flow U, B is the magnetic 
field normalised by v'47rp, p = 1 is the density, p is the 
pressure determined by requiring V • n = 0, and v and r] 
are the kinematic viscosity and magnetic diffusivity coeffi- 
cients, respectively. 

Using this local formulation allows us to study differ- 
ent instances of differentially rotating background flows by 
changing the values of the shearing and rotation rates S and 
ft. In this paper, we consider two important cases: (a) shear 
with no rotation, S y^ 0, ft = (see Sect. O and (b) Ke- 
plerian rotation, Q = —2S/3 (see Sect. |4]). The nonrotat- 
ing ca se has been discussed by us previously (I Yousef et al. 



2008h and is relevant to nonrotating systems as well as 



systems where one expects the shear to be stronger than 
the rotation (one good astrophysical example is irregular 
galaxies, where large-scale magnetic fie lds are found in the 
absence of strong overall rotation Ij see Chvzv et al.l l200d 



I2OO3I: iGaensleretallbOOSl: iKeplev et al.ll2007k The case of 
Keplerian rotation corresponds to the local shearing- sheet 
model of thin accretion disks, for which the angular velocity 
decreases radially ex R~^^'^. It is commonly used to study 
turbulence and transport in accretion disks, in particular the 
MHD turbulence drive n by the magneto-rotational insta- 
bilitv (MRP (see, e.g.. rBalbusll2003riFromang et"aD 12007 : 
Lesur & Longarettill2007[ and references therein). 



We solve Eqs. ([B and Q in a shearing periodic compu- 
tational domain by means of a L agrangian spectral method 
dLithwickl l2007[ lOgilvid [19981) . In our model, we use a 
random body force / that is nonhelical satisfying / • 
(^ >< /) = O7 (5-correlated in time (white noise), and has 
a characteristic length scale // = 27r//c/. Choosing / to be 
^-correlated in time makes it possible to impose the aver- 
age injected power e = {u ■ f) precisely. Here and in what 
follows, angular brackets denote volume averaging, while 
overbars represent time averaging. The characteristic outer 
scale length of the turbulence excited by our random body 



We thank D. Sokoloff for pointing this out to us. 
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Table 1 Index of runs. 
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Q 


L. 


Resolution 


v = 


v 


7 


Ib 




Run 


B</B< ^ 


SI 


-2 





8 


32^ X 256 


1 X 


10-2 


0.0161 


3.7 


11.5 


S2 


-2 





16 


32^ X 512 


1 X 


10-2 


0.021 


3.8 


11.1 


S3 


-1 





8 


32^ X 256 


1 X 


10-2 


0.0030 


4.6 


10.0 


S4 


-1 





16 


32^ X 512 


1 X 


10-2 


0.0124 


5.4 


10.2 


S5 


-1 





32 


32^ X 1024 


1 X 


10-2 


0.0092 


5.2 


9.9 


S6 


-1 





64 


32^ X 2048 


1 X 


10-2 


0.0121 


5.1 


9.5 


S7 


-1/2 





16 


32^ X 512 


1 X 


10-2 


0.0040 


6.8 


9.1 


S8 


-1/2 





32 


32^ X 1024 


1 X 


10-2 


0.0058 


7.1 


9.0 


S9 


-1/2 





64 


32^ X 2048 


1 X 


10-2 


0.0055 


7.3 


9.0 


SIO 


-1/4 





64 


32^ X 2048 


1 X 


10-2 


0.0025 


9.7 


8.3 


Sll 


-1/4 





128 


32^ X 4096 


1 X 


10-2 


0.0025 


9.9 


8.2 


S12 


-1/8 





64 


32^ X 2048 


1 X 


10-2 


0.00094 


13.1 


8.1 


S13 


-1/8 





128 


32^ X 4096 


1 X 


10-2 


0.00092 


13.5 


8.1 


Kl 


-1/2 


1/3 


32 


32^ X 1024 


1 X 


10-2 


0.0067 


7.0 


8.9 


K2 


-1 


2/3 


16 


32^ X 512 


1 X 


10-2 


0.0174 


5.0 


9.7 


K3 


_2 


4/3 


8 


32^ X 256 


1 X 


10-2 


0.038 


3.5 


10.7 


K4 


-4 


8/3 


8 


32^ X 256 


1 X 


10-2 


0.086 


2.6 


12.7 


K5 


-8 


16/3 


8 


32^ X 256 


1 X 


10-2 


0.118 


2.2 


20 


K6 


-16 


32/3 


8 


32^ X 256 


1 X 


10-2 


0.147 


2.1 


44 


K7 


-32 


64/3 


8 


32^ X 256 


1 X 


10-2 


0.078 


2.3 


137 


K8 


-64 


128/3 


8 


32^ X 256 


1 X 


10-2 


0.036 


2.9 


520 


FDl 


-1 


2/3 


8 


64^ X 1024 


1 X 


10-3 


0.32^ 


5.8^^ 


16.4^ 


FD2 








8 


64^ X 1024 


1 X 


10-3 


0.29^ 


- 


- 














values differ from the 


corresp 






is actual! 








)rted \ 




^ This 


y [{B<^)./{B<'- 


^>^]i/2.Therepc 


onding numbers 



(By)^ I {Bx)xy (fields averaged over x and y) instead of [{By )z/ {Bx )zY^'^ (^s erroneously claimed there). 



y I xy 

^ Calculated for the kinematic regime of the fluctuation dynamo (where Brms /^r 
and Brms for runs FD1-FD2. 

^ Calculated for period after saturation of the fluctuation dynamo (t > 200). 



; < 10 ^). Refer to Fig. [9] (left panel) for time evolution of i 



force is ~ // and for a given viscosity, the corresponding 
outer scale velocity ~ i^rms = {u^Y^'^ is then determined 
by e; for sufficiently small viscosities, dimensionally, we 
must have Urms ^ {^hV^^- 

In astrophysical systems, where the source of the tur- 
bulence is often related to the large-scale shear itself, the 
most relevant case tends to be iS'r ~ 1, where r ~ If/urms 
is the turnover time of the turbulent motions at the outer 
scale. In our model, we have control over the turbulence via 
the parameters of the body force described above. This al- 
lows us to consider a range of values of S. Such a study is 
useful because 5'r ^ 1 is only an order-of-magnitude rela- 
tion (so it is not obvious which precise value of S would be 
the most "realistic") and also because any analytic theory 
of the shear dynamo would have to predict how the effect 
depends on the value of S (usually it is only the asymptotic 
case Sr <^ 1 that can be treated an alytically, as, e.g., in 
iRogachevskii & Kleeorinll2003ll2004l) . 

As we shall see, a significant separation in scales be- 
tween If and the dimensions of the computational domain 
is necessary to get magnetic field growth — this is partic- 
ularly true for weak shears. Also very long run times (at 
least hundreds and, for weak shears, thousands of turnover 
times) are needed to get reliable statistics. This is compu- 



tationally expensive and we were only able to satisfy both 
requirements within available computational resources by 
restricting our study to computational domains elongated in 
the z direction. Thus, in all the runs presented here, we have 
Lx = Ly = 1, the forcing scale is // = 1/3 and Lz is taken 
to be between 8 and 128 (depending on the value of S), so 
Lz ^ //. The minimum value of L^ needed for a given 
simulation to produce results independent of the box size in 
the z direction depends on the shearing rate (this issue is 
discussed in Sect. l3.3l) . Setting Lx = Ly = 1 together with 
e = 1 defines the code units of length and time. 

A list of all runs is given in Table [U The S (sheared, 
nonrotating) and K (Keplerian) runs have v = r] = 10"^. 
The S runs have u^^s = 1.0 ± 0.3 while the K runs have 
^rms = .5 ± 0.06 (larger velocities in the former case are 
due to the excitation of large-scale velocity structures in 
the absence of rotation; see Sect. 13.41) . This corresponds to 
Reynolds numbers Re = Urms/kfi' ~ 5 and 3 for the S 
and K runs, respectively. The magnetic Reynolds number 
Rm = Urms/kff] = Re for all runs, as the study of the de- 
pendence of the shear dynamo properties on the magnetic 
Prandtl number Pm = Rm/Re lies outside the scope of 
this paper. As the Reynolds numbers are low, only relatively 
low resolution is needed and we find it sufficient to have 32 



www.an-journal.org 



© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



740 



T. A. Yousef et al.: Numerical experiments on dynamo action in sheared and rotating turbulence 




Fig. 1 (online colour at: www.an-journal.org) The upper 
panel shows the evolution of 5rms and Uj-ms with time for 
the run K2 (S = —1,1] = 2/3). The lower panel shows 
the normalised kinetic helicity as a function of time t for 
the same run. Fluctuations in kinetic helicity only occur on 
turbulent time scales and have zero time average. This is 
true for runs both with and without rotation. 



collocation points per code unit of length. Strictly speak- 
ing, this is not fully developed turbulence as there is little, if 
any, separation between the energy-containing (outer) and 
viscous scales. However, this is not a serious limitation be- 
cause one expects the mean field dynamo mechanism and 
the properties of the resulting mean field to depend only on 
the outer scale of the turbulence (rather than on the structure 
of the inertial-range) and, therefore, to become asymptoti- 
cally independent of Re for some critical Re > 1 While 
one cannot, of course, always take this kind of intuition for 
granted, one expects it to work at least for a turbulence 
driven with fixed input power e as in our simulations. It 
does, indeed, appear to work: by rerunning a selection of 
the runs with lower viscosity, we have confirmed that the 
results reported below are independent of Re. 



All simulations start with a dynamically insignifi- 
cant random seed magnetic field whose volume aver- 
age is (and remains) zero. The initial magnetic energy is 
(52)^10-20 2/2^^ for runs S1-S13 and (^^^-IQ-^^ ^^r 
for runs K1-K8 and FD1-FD2. The shear dynamo mecha- 
nism generates a magnetic field with large-scale structure 
in the z direction. The field grows exponentially until it 
reaches a saturated state, in which the value of {B'^) be- 
comes comparable to {u'^). In this paper, we only analyse 
the kinematic regime of this dynamo, in which the back re- 
action of the magnetic field onto the velocity field through 
the Lorentz force can be neglected, and so only consider 
data for times before {B'^) reaches 10~'^{u'^). All numbers 
reported in Table \T\ are for this kinematic stage, except for 
the run FDl (see Sect.©. 



Before we close this section, several important remarks 
are in order, concerning the physical effects that are not 
present in our simulations. 

First, the magnetic Reynolds numbers Rm for all runs 
except FDl and FD2 are deliberately chosen to lie below 
the critical threshold value Rm^fd at which the fluctuation 
dynamo sets in. This ensures that the magnetic-field growth 
observed in the S and K runs is a result of a "pure" mean- 
field dynamo, unpolluted by the fluctuation dynamo and that 
the operation of the former is not predicated on the existence 
of the latter. 

Second, the magnetorotational instability, MRI turbu- 
lence or MRI dynamo are not present in the Keplerian runs 
K and FDl. This is because the magnetic field is extremely 
weak (in the K runs) and the viscosity is large. In addition, 
we know that the MRI dynamo and turbulence have so far 
only been found for values of Pm systematically larger than 
unity, in high-resolution simulations of Keplerian shearing 
boxes elongated in th e y direction (IFromang et al.l l2007l: 
Lesur & Ogiivi3l2008h . 

Third, the turbulence in all our simulations is statisti- 
cally nonhelical, even for the runs with rotation. The ki- 
netic helicity Hk = (n • (V x n)) has a zero time av- 
erage, Hk = 0, and only fluctuates on the turbulent time 
scales (~ r ~ If/urms)- This is a consequence of / being 
nonhelical, / • ( V x /) = 0, and (5-correlated in time. The 
latter guarantees that the time average of helicity production 
is ( V X n) • /+ n • ( V X /) = 0. Figure [T] shows the time 
evolution of Hk for a typical run. This excludes the possi- 
bility that a effect, at least in its standard and simplest form, 
is present in our simulations. 



3 Shear dynamo without rotation 

In this section, we document the numerical results obtained 
for 5* 7^ and 17 = (turbulence + shear, no rotation). 
This regime is covered by the S runs (Table [T]). In all of 
them, we find that magnetic field grows exponentially. Let 
us first describe the spatial structure of this growing field, 
then the dependence of its properties on the shearing rate S 
and their convergence with respect to the vertical size of the 
computational box. 



rms 3.1 Structure and evolution of the growing field 

Figure [2] shows a cross section of the different components 
of B and u from a representative run with no rotation and 
weak shear. One can clearly distinguish between fluctua- 
tions with scales comparable to the forcing scale (which 
is, on the average, // = 1/3 of the horizontal box size 
Lx = Ly) and a large-scale "mean" field. The small-scale 
fluctuations cannot result from the fluctuation dynamo be- 
cause Rm < Rmc,fd. Instead, they stem from the tangling 
of the mean field by the turbulence - the energy in these 
fluctuations is comparable to the energy of the mean field 
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Fig. 2 (online colour at: www.an-journal.org) Snapshots of Bx, By, Bz, Ux, Uy and Uz in the (?/, z) plane for run S5 
{S = —1,Q = 0) at t = 300. Both velocity and magnetic fields have structure in the z direction with characteristic 
scale far larger than the horizontal box size Lx = Ly = 1, in addition to the turbulent fluctuations that have characteristic 
length scales smaller than 1 (see Sect. [331 regarding the velocity structure). The panels underneath the snapshots show the 
instantaneous profiles of the field components (for a given x = xq and y = yo) vs. z (thin lines) and their large- scales 
means Fourier-filtered according to Eq. Q (thick lines). 



(for larger Rm, this ene rgy will be even larger; see, e.g., 
ISchekochi hin et aLOOOVl and references therein). 



iqo f 



The spatial structure of the magnetic field is illustrated 
by Fig.O which shows the time-averaged normalised spec- 
trum of the growing field with respect to kz , 



A^(M = ^ E \Bit.C,ky,k.W 



(3) 



l^x il^v 



Note that since the average forcing wave number is kf = 
{klf + klf + kljY'^ = 27r/lf, where Ij = 1/3, the 
mean square forcing wave number associated with one spa- 
tial dimension is kzj = (A:?/3)^/^ = 27r>/3. The peak of 
the spectrum at kz <C kzj corresponds to the mean field. 
The fluctuations at smaller scales are also evident, although 
at the low value of Rm that we use they are not very large. 

In order to systematically isolate the mean field from 
the small-scale fluctuations, we filter out the highest Fourier 




kz/27r 

Fig. 3 (online colour at: www.an-journal.org) Time aver- 
aged normalised one-dimensional spectra of the magnetic 
energy (Eq. 3) for runs S3-S6 (S = l,ft = 0, vertical box 
sizes Lz = 8, 16, 32, 64). As Lz is increased, the peak of 
the spectrum settles at a wave number independent of L^ 
(see Sect. [33]). 
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modes: 



/e^/27r<l 



ikzZ 



(4) 



This is our definition of the mean field. The Fourier filtering 
averages out all variation in the x and y directions and re- 
moves the scales in the z direction smaller than L^ = Ly = 
1 = 3//. Note that B^ is strictly zero because V • ^ = 0: 
indeed, d^B< = -5^5< - dyB< = as the filtered field 
has no variation in x and y. Of the two nonzero components 
of the growing mean field, B^ component is the larger: 
this can be seen in Fig. [21 where we plot the instantaneous 
profiles of the field components and their Fourier- filtered 
means. This relative predominance of By is not surprising 
given the form of Eq. Q where in addition to the turbulent 
stretching terms, there is a linear term BxSy whereby the 
shear systematically produces By from Bx (azimuthal field 
from radial field in axisymmetric systems; this is known as 
the l]-effect in dynamo theory). 

As well as having a large-scale spatial coherence, the 
mean field is also coherent over long times. While it does 
not seem to be an eigenmode in the usual sense of having 
a shape that is exactly constant in time, the characteristic 
times over which it changes are much longer than both the 
turbulent turnover time r and the shearing time S~^. The 
time evolution of the mean field is illustrated by Fig. IH (the 
left panels for nonrotating runs, the right panels for rotating 
runs; see Sect.©. Here we plot the evolution of the profile of 
the y component of the magnetic field averaged over x and 
y, (By)^ {Zjt). In all cases, we see a well-defined evolving 
spatial structure: patches of stronger or weaker field slowly 
drift along z, occasionally peter out and reemerge. 

3.2 Dependence on the shearing rate 

Dependence of the properties of the growing field on the 
shearing rate S is the first obvious parameter scan to carry 
out. It provides a quantitative test of past mean-field theories 
and possibly a valuable hint as to how one might construct 
future ones. 

The dependence of the growth rate 7 of Brms = 
{B^y^^ on S is plotted in the left panel of Fig. [3 The 
growth rates for the nonrotating runs appear to fit approx- 
imately the dependence 7 oc 6* for 5* between 1/8 and 2. 
Note, however, that the actual numerical values of 7 are 
substantially smaller than S (see Table [T]). We do not cur- 
rently know what determines the proportionality coefficient 
(we have checked that it is not the value of Re: the growth 
rates presented here do not change as the Reynolds number 
is increased). 

In Tabled we also fist [/ B< '^dz/ J B<'^dz]^. Its value 
is approximately 10 irrespective of S , although it shows a 
slight tendency to increase with S. This is consistent with 
J (X S because from the y component of the induction equa- 
tion Eq. ([2|, we can estimate dfB^ ^ jB^ ^ SB^, so 
B</B< -5'/7- const. 
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Fig. 6 (online colour at: www.an-journal.org) Bottom 
panel: time evolution of 5rms for a set of runs with S = \ 
and varying vertical box size L^ = 8, 16, 32, 64 (runs S3- 
S6). When L^ is sufficiently large (here 16) the growth rate 
7 becomes independent of L^. Top panel', evolution of ^/rms 
for the same runs. The velocity fields develops large long- 
lived fluctuations (see Sect. 13.41) . 



From Fig.m we already knew that the scale of the mean 
field decreases with shear. In order to study this dependence 
quantitatively, we define the characteristic length scale of 
the mean field as follows 



1 



jdz{dB</dzy 

JdzBf 



1/2 



(5) 



The values of Ib for all runs are listed in Table [T] and plotted 
in the right panel of Fig. [5] as a function of the shearing 
rate S. We see that for the values of S that we have studied, 
the characteristic length scale Ib approx imately matches the 
scaling 1/y/S. It is possible to argue (I Yousef et al.l l2008k 
that for a generic model form of mean-field equations, this 
scaling would indeed be consistent with the linear scaling 
of the growth rate, j oc S. 

These scalings do not agree with the mean field the- 
ory in its current state. Namely, iRogachevskii & KleeorinI 
(|2003|) . that proposed the existence of the shear dynamo us- 
ing a mean-field theory with the r-approximation closure to 
incorporate the effect of turbulence, predicted that for the 
fastest-growing mode 7 ex S*^ and Ib cx 1/S for S'r <C 1. 
These scalings do not seem to match the numerical results. 
However, we cannot exclude the possibility that the values 
of S considered by us are, in fact not asymptotic in the small 
parameter Sr and that Rogachevskii and Kleeorin's scalings 
might hold for even weaker values of shear or for numeri- 
cal domains that allow more scale separation between the 
turbulence and the horizontal box size (Lx and/or Ly). 

3.3 Independence of the computational domain size 

It is important to ensure that the results reported above, par- 
ticularly those of Fig. [Jl are independent of the size of the 
computational domain. Since the growing mode has large- 
scale structure in the z direction, we do this by performing. 
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Fig. 4 (online colour at: www.an-journal.org) Left panel: the y component of the magnetic field averaged over x and y, 
{By)^ (2:, t), for three runs with linear shear and S = 1/2,1,2 and 17 = (runs S8, S4 and SI, from top to bottom). 
The characteristic length scale Ib (manifested by variations along the vertical axis of the figures) decreases with increasing 
shearing rate S. The magnetic structures are correlated over times (variations along the horizontal axis of the figures) that 
are very long compared to the turnover time of the turbulence r ~ 1/3. Right panel, similar figures for the case with 
Keplerian rotation (runs Kl, K2, K3) discussed in Sect. HI 
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Fig. 5 (online colour at: www.an-journal.org) Left panel, the growth rate 7 of B^ms vs. the shearing rate S. The numerical 
values of 7 are given in Table [B For each value of S, only the run with the largest value of Lz is shown in the figure to 
avoid clutter. The plotted values correspond to the nonrotating runs S2, S6, S9, SI 1, S 13 (red circles; see Sect. l3.2l) and the 
Keplerian runs K1-K8 (green squares; see Sect. 131). Right panel, the characteristic scale Ib of the mean field, as defined by 
Eq. (|5]) vs. the shearing rate for the same runs (see Table[T]for the numerical values of Ib)- Note that only small shear values 
are accessible in the nonrotating case because larger shears encounter nonlinear instability in that regime (see Sect. 13.41) . In 
the Keplerian regime, these nonlinear instabilities are suppressed, which makes it possible to explore regimes with larger 
shear (see Sect.©. 



for each value of S, a set of runs in boxes with successively 
larger values of L^ until we can ascertain that the proper- 
ties of the growing mode are independent of Lz. To illus- 
trate this procedure, the bottom panel of Fig. [6l shows the 
time evolution of B^ms for runs with 5 = 1: exponential 
growth with a growth rate independent of L^ emerges as L^ 
is increased. Another illustration of the convergence with 
respect to L^ is Fig. [3l where normalised magnetic-energy 
spectra are plotted for the same runs as in Fig.O We see that 



the peak of the spectrum (corresponding to the mean field) 
becomes independent of Lz at sufficiently large Lz. 

We have performed this type of study for a number of 
values of S and in all cases confirmed convergence with re- 
spect to Lz. This is documented in Table [TJ Naturally, as 
the scale of the mean field decreases with increasing shear- 
ing rate S (Sect. [T2l) , the minimum Lz required for conver- 
gence also decreases. Thus, we find that Lz = 8, 16, 32, 64 
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are sufficiently large for shearing rates S = 2, 1,1/2, 1/4 
and 1/8, respectively. 

The convergence with respect to L^ ensures that the 
properties of the mean field and their dependence on the 
shearing rate S are physical and are not numerical artifacts 
related to the choice of the computational domain size in z. 
It is also important to check whether the results depend on 
the dimensions of the computational domain in the x and 
y direction. Such study, which is more computationally de- 
manding, has been left outside the scope of this paper and 
will be undertaken elsewhere. 

3.4 Vorticity dynamo 

We finally mention another noteworthy result in the nonro- 
tating regime. We observed that the velocity field also de- 
velops energetic fluctuations with correlation times much 
greater than r and S~^. This property is illustrated in the top 
panel of Fig. [6l Since homogeneous linear shear is known 
to become nonlinearly unstable (see Pumirl ll996[ and ref- 
erences therein) at large enough Re, one might imagine 
that the observed velocity fluctuations result from finite- 
amplitude destabilisation of shear by forcing- scale motions. 
Indeed, for high enough values of S and/or L^ we experi- 
enced that our background shear became vigorously unsta- 
ble and developed into a flow characterised by Urms ^ 1 
and fluctuations on all scales. However, this is not the case 
for any of the runs presented here. There is a key dif- 
ference between these runs and the truly unstable cases. 
Finite-amplitude inst abilities are known to feed on shear 
(ICasciola et al.ll2003|) and are characterised by large (uxUy) 
correlations associated with turbulent momentum transport. 
This contrasts with our results: as Fig.|7]shows, there seems 
to be almost no correlation between the large velocity fluc- 
tuations and the mean power injected into the flow by 
the shear —S{uxUy) in our runs. The ratio of this shear 
power to the forcing power e stays much smaller than unity, 
—S{uxUy)/e <c; 1, implying that the velocity fluctuations 
feed on the forcing and not on the shear. 

These observations suggest that the finite- amplitude in- 
stabilities associated with the shear are not the cause of the 
observed velocity fluctuations. Given the formal similarities 
between the induction and vorticity equations, one might ar- 
gue that a ' Vorticity dynamo" similar to the shear dy namo 
is triggered (lElperin. Kleeorin & Rogachevskiill2003|) . This 
is supported by the fact that the large-scale velocity struc- 
ture forms mainly in Uy (see Fig. [2]), which corresponds to 
the large-scale vorticity ujx = —dzUy, the quantity one ex- 
pects to be analogous to By. The difference between the 
vorticity and the magnetic-field dynamos is that the "seed" 
field for the vorticity dynamo results from the forcing and 
is, therefore, never dynamically insignificant. Thus it is not 
possible to observe a long kinematic growth phase for this 
vorticity dynamo. One might conjecture that the large ve- 
locity fluctuations represent the nonlinear regime of such a 
dynamo. We note, however, that the status of the vorticity 
dynamo remains uncertain because, as seen, e.g., in Fig.[6l 
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\ / 




Fig. 7 (online colour at: www.an-journal.org) Top panel: 
time evolution of u^rns for the runs S4 (5 = — 1, ^^ = 0; 
dashed line) and K2 (5 = — 1, 1^ = 2/3; solid line). Large- 
amplitude long-lived in Urms can be seen for run S4, but not 
for run K2. Bottom panel, evolution of the power injected 
into the flow by shear —S{uyUx) for the same runs. This 
is to be contrasted with the power injected by the forcing, 
e = 1. 



both the time duration of the large velocity fluctuations and 
their amplitude depend strongly on L^. Thus, unlike for the 
magnetic field, there is no convergence here with respect to 
the domain size and we must remain cautious in interpreting 
these results. 

Finally, we observe that the emergence of the large- 
amplitude, long-lived, large-scale fluctuations in the veloc- 
ity field does not appear to be strongly correlated with the 
operation of the shear dynamo: compare, e.g., the time evo- 
lution of iirms and 5rms shown in Fig.[6l As we are about to 
see, the absence of these velocity fluctuations in the simu- 
lations with Keplerian rotation (see Fig. [7]) does not change 
any of the basic properties of the shear dynamo, so we feel 
safe in ruling out the possibility that the large-scale velocity 
fluctuations are a required ingredient in the shear dynamo. 

4 Shear dynamo in rotating systems 

Since many astrophysical plasmas (stars, disks, galaxies) 
are differentially rotating, it is important to understand 
whether, and how, the shear dynamo is affected by im- 
posing a uniform rotation of the system. We focus on the 
case of Keplerian rotation introduced in Sect. O because 
it is relevant for accretion-disk dynamos and has b een ex- 
tensively studied numerical ly in this context (e.g., iBalbus 
20031: iFromang et al.l l2007l and references therein). This 
regime is thought to be nonlinearly stable from the hydro- 
dynamic point of view ( Ji et al.l 20061: Lesur & Longaretti 

:ossul I2007I) . 



20051: iRincon. Ogilvie & Cossul I2007h . meaning that the 
nonlinear instabilities characteristic of nonrotating shear 
flows are absent. They are, indeed, absent in our simula- 
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tions. Furthermore, as illustrated by Fig. [71 even the large 
velocity fluctuations present in the time series of Urms for 
the nonrotating runs are suppressed in the case of Keplerian 
rotation. The suppression of these fluctuations, which were 
tentatively interpreted in terms of a "vorticity dynamo" in 
Sect. 13.41 is likely to be another consequence of the stabilis- 
ing effect of Keplerian rotation on hydrodynamic motions. 



As we have mentioned in Sect. 13.41 the nonlinear desta- 
bilisation of the shear for too large values of S and Lz im- 
posed limitations on the range of values of S for which shear 
dynamo could be studied in the nonrotating case. The ab- 
sence of such destabilisation in the simulations with Keple- 
rian rotation allows us to extend the study the shear dynamo 
to much higher values of S than in the nonrotating case (see 
the K runs in Table [T]). 

The main conclusion of this study is that Keplerian ro- 
tation does not seem to alter the properties of the shear 
dynamo in any significant way. All the results concern- 
ing the growing magnetic field reported in Sect. [51 both 
qualitative and quantitative, continue to hold in the rotat- 
ing case. The right panel of Fig. |4l shows the evolution of 
the large-scale field with time for three runs with rotation. It 
is hard to see any qualitative difference between the struc- 
ture and evolution of the growing field in the rotating case 
and its structure in the nonrotating case shown in the left 
panel of the same figure. The similarity between the two 
cases is further demonstrated in Fig. \5\ which shows that 
for low values of S, both the growth rates 7 and the char- 
acteristic length scales Ib seem to follow the same scaling 
laws in the rotating and nonrotating cases. The numerical 
values of 7 with and without rotation are quite close, al- 
though the growth rates in the presence of rotation appear 
to be systematically slightly higher. This small difference 
may result from the changes in the structure of the tur- 
bulence caused by rotation or it could be due to an addi- 
tional contribution to the generation of magnetic field from 
other mean-field dynamo mechanisms assoc iated with ro- 
tation: e.g., the Radler, or ft x J, ef fect dRadled fl969l: 
iRadler. Kleeorin & Rogachevskiill2003b . It is worth reem- 
phasising in this context that the turbulent motions driven by 
our forcing remain nonhelical on average even in the pres- 
ence of Keplerian rotation (see Fig. [T] and the discussion at 
end of Sect. [J]), so that the rotating shear dynamo does not 
involve an a effect in the usual sense. 

Figure \5\ shows that, as S is increased beyond the val- 
ues studied for ^^ = 0, the scalings of the growth rate and 
the mean-field scale, 7 oc 6* and Ib cx l/V^, continue un- 
til 7 reaches its maximum and Ib its minimum at 6* ^ 10. 
At even larger shears, the growth rate starts to decrease, as 
it intuitively should. Indeed, the dynamo must surely disap- 
pear in the limit of infinite shear. More precisely, in the limit 
when the shearing rate associated with the imposed shear is 
much larger than the turbulent rate of strain, S :^ \ Vn | , the 
nonlinear terms (the terms involving turbulent velocities) in 
the induction equation Eq. (O become negligible — without 




kz/27r 

Fig. 8 (online colour at: www.an-journal.org) Time aver- 
aged normalised one-dimensional spectra of magnetic en- 
ergy (Eq. 3) for runs K1-K5 (S = -\, -1, -2, -4, -8; 
Vl = — 2S'/3). As S decreases, the peak of the spectrum, 
corresponding to the mean field, moves towards lower wave 
numbers, corresponding to larger scales. 



these nonlinearities, we cannot have an exponentially grow- 
ing magnetic field. 

This argument suggests that the qualitative change asso- 
ciated with the peak of the growth rate in Fig.[5]corresponds 
to the transition from the weak- shear regime {St <C 1) to 
the strong- shear regime {St ^V). Quantitatively, this tran- 
sition evidently occurs at 5* ~ 10. Once the imposed shear S 
becomes stronger than the turbulent rate of strain at the forc- 



ing scale r 



~ Ur 



s/lf, smaller turbulent scales should 



come into play. In order to understand what happens when 
the shear is strong, one must analyse the changes that such a 
shear causes to the inertial-range turbulence, determine the 
transition scale at which the turbulent rate of strain catches 
up with the shear, etc. Since the numerical simulations pre- 
sented here, unlike real turbulent systems, have no extended 
inertial range, their relevance in the strong-shear regime is 
questionable and we shall not discuss this regime here. 

Finally, for completeness, in Fig. [8] we show the one- 
dimensional magnetic-energy spectra, defined by Eq. ([3]), 
for Keplerian runs with several values of shear. We see that, 
as S is decreased, the peak of the spectrum, correspond- 
ing to the mean field, moves towards lower wave numbers 
(corresponding to larger scales; cf. Fig. [5]). The rest of the 
spectrum represents magnetic fluctuations produced by the 
turbulent tangling of the mean field. Here, like for the non- 
rota ting runs in Sect. [3 the ta ngling, or magnetic induction 
(see lSchekochihin et al.l2007[ and references therein), is the 
only mechanism for producing small-scale fields because 
Rm is chosen below the threshold for the fluctuation dy- 
namo. 

5 Shear dynamo and fluctuation dynamo 

So far, we have only discussed simulations with Rm be- 
low the critical value for the onset of the fluctuation dy- 
namo Rmc fd- This choice was made to ensure that the 
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Fig. 9 (online colour at: www.an-journal.org) Left panel: time evolution of Brms for runs FDl (5 = — 1 and Q = 2/3) 
and FD2 (S = ^ = 0). The inset shows the same evolution on a linear scale. Right panel: The y component of the magnetic 
field averaged over x and y, {By) (z, t), for runs FDl (left) and FD2 (right). The run FDl develops large-scale structure 
with long correlation time, while the run FD2 only develops small short-lived fluctuations. 




rate of strain associated with scales much smaller than the 
outer scale, so 7 :^ r~^, where r is the outer- scale tur- 



kj27T 

Fig. 10 (online colour at: www.an-journal.org) Time aver- 
aged normalised one-dimensional magnetic energy spectra 
(Eq. 3) for FDl (S = -lM = 2/3) andFD2 (S = Q = 0). 
The time average is taken over t = 200 to 500. Unlike for 
FD2, the spectrum for FDl peaks at low wave numbers (cor- 
responding to large scales). This indicates the existence of 
large-scale magnetic structure. 



growth in magnetic energy is solely due to the shear dy- 
namo, allowing us to confirm the existence of this effect, 
its independence of the fluctuation dynamo and to study 
it in a "pure" isolated form. However, astrophysical ob- 
jects virtually always have large Rm, which are expected 
to be far above the fluctuation dynamo threshold. Thus, 
any mean-field dynamo scheme aspiring to explain large- 
scale astrophysical fields has to work simultaneously with 
the fluctuation dynamo. Furthermore, the fluctuation dy- 
namo is likely to be much faster than any mean-field dy- 
namo because its growth rate is determined by the turbulent 



bulent turnover time (see discussion in ISchekochihin et al.l 
20071 and references therein). In contrast, mean-field dy- 
namos are typically slower than the outer-scale rate of 
strain, so for them, 7 <C r~^. Thus, a working mech- 
anism for mean-field generation must operate against the 
background not of pure hydrodynamic turbulence but rather 
of magnetohydrodynamic turbulence that emerges as the 
saturated state of the fluctuation dynamo (the approach 
advocated, e.g.. bvlRadler. Kleeorin & Rogachevskiill2003l: 



Rogachevskii & KleeorinI 120031 l2004l) . It cannot be taken 
for granted that such a magnetohydrodynamic turbulence in 
conjunction with shear will still be germane to the gener- 
ation of large-scale magnetic fields. It is, therefore, imper- 
ative to confirm the existence of the shear dynamo effect 
for Rm > Rmc,fd in order to be able to assert its poten- 
tial relevance for astrophysical situations. Here we show (as 
far as we know, for the first time) that a combination of lin- 
ear shear and a nonhelically forced homogeneous magneto- 
hydrodynamic turbulence (saturated state of fluctuation dy- 
namo) can support dynamically significant large-scale mag- 
netic fields. 

For this preliminary study, we performed two runs, FDl 
and FD2. For FDl, S = -l, 9. = 2/3; for FD2, S = 
ft = 0. Both runs have u = j] = 1 x 10~^, iXrms ^ 1 
and thus Re = Rm ^ 50. This value is above the thresh- 
old for the fluctuation dynamo. In the run FD2, which is a 
pure homogeneous turbulence run, that is the only possible 
field-amplification mechanism. In contrast, run FDl has the 
prerequisites for both the shear and fluctuation dynamos. 

The time evolution of the magnetic field for these two 
runs is shown in Fig. [51 Both runs develop in a sim- 
ilar manner initially: the magnetic field is amplified by 
the fluctuation dynamo and 5rms increases by 12 orders 



© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2008) 



747 



of magnitude to a dynamically strong saturated level af- 
ter a time t ~ 100 ~ lOOr. This growth is much 
faster than the growth rates associated with the shear dy- 
namo in analogous runs without the fluctuation dynamo. 
The growth rate is 7 c^ 0.3, which is consistent with 
the expectation that the growth rate of the fluctuation 
dynamo should be comparable to the turbulent rate of 
strain, or stretching rate — for Rm > Re, this is usu- 
ally the rate of strain associated with the smallest turbu- 
lence scales (see, e .g., 'Bran denburg & SubramanianI 12005 : 
ISchekochihin et al.ll2004. .20071 and references therein), but 
in our case the scale separation between the outer and the 
viscous scale is not large, so the fluctuation-dynamo growth 
rate should be roughly comparable to r~^. 



As evident in Fig.[9l the two runs differ significantly af- 
ter the saturation of the fluctuation dynamo. The magnetic 
field for the zero- shear run FD2 becomes statistically sta- 
tionary and stays at small scales, while in the sheared and 
rotating run FDl it develops large-scale structure. The struc- 
tural difference between the runs with and without shear 
is displayed in the right panel of Fig. O whereas the zero- 
shear run only has small scale, short-lived turbulent fluctua- 
tions, the run with an imposed shear develops large-scale, 
long-lived structures, which are in many ways similar to 
the shear-dynamo-generated mean field in the runs with- 
out the fluctuation dynamo. The existence of large-scale 
structures in the case with shear can also be seen from 
(time-averaged normalised one-dimensional) magnetic en- 
ergy spectra shown in Fig.[TOl The spectrum for FDl (S = 
-l,n = 2/3), unlike that for FD2 (S = n = 0), has a 
prominent peak at a low wave numbers corresponding to 
a large-scale mean field. Note that 5rms in FDl appears to 
develop long-lived, large- amplitude, large-scale fluctuations 
that are reminiscent of those observed in Urms for the vor- 
ticity dynamo reported in Sect. 13.41 (see the inset in the left 
panel of Fig. [9]) . 



Thus, we have found that the fluctuation dynamo and 
the shear dynamo can coexist, giving rise both to small- 
and large-scale magnetic fields. It is clear that a system- 
atic parameter study at high Rm and Re is called for. It 
requires much larger computational resources than the rel- 
atively cheap simulations reported in this paper. Since the 
seed field for generating large-scale magnetic fields is now 
provided by the already saturated small-scale fields result- 
ing from the fluctuation dynamo, it is never really well 
posed under these circumstanc es to consider a kinemati c 
regime of the shear dynamo (cf . ICattaneo & Hughesll2008|) . 
Thus, we have to confront the question of how the large- 
scale fields saturate and what their structure is in a fully 
nonlinear dynamical regime. We have left these questions 
outside the scope of this paper. Further investigations along 
these lines will be reported elsewhere. 



6 Conclusion 

Can nonhelical turbulence in combination with large-scale 
velocity shear act as a mean-field dynamo and generate 
magnetic fields with length scales much larger than the outer 
scale of the turbulen ce? This has been the subject of consid- 
erable r ecent deb ate ^Brandenburg 2005; Brandenburg et aL 




iFedotovl 2003; Fedotov, Bashkirtseva & Rvashkg 
Kleeorin & Rogachevskiil I2OO8I: IProctod l2007l: 



Radler & Rheinhardt' '2OO' 



20061: JRogachevsk i i & Kleep rii 



eeorini 



Rudiger & Kichatinovl I2OQ6I: ISilant'ev 
1999a||bl I2OO2I; IVishniac & Brandenburg' 




To our 

knowledge, the paper by Yousef et al. (2008) and this paper 
present the first set of dedicated numerical experiments 
that demonstrates that such a generation mechanism is 
feasible. However, in retrospect, one might conjecture 
that the shear-dynamo might have already been seen in 
several earlier numerical studies that combined large-scale 
flows, and consequently large-scale shear, with nonhelical 
forcing at smaller (or, in some cases, similar) scales and 
reported generation of magnetic fields at scales larger than 
the forcing scale (Brandenburg 200 5; Mininni et al. 2005|; 
Pontv et al.l[2005l: IShapovalovll2006h . 

We have carried out a suite of numerical experiments 
on the shear dynamo effect in vertically elongated shear- 
ing boxes and for magnetic Reynolds numbers subcritical 
with respect to the fluctuation dynamo. For the values of 
the imposed shear S between 1/8 and 8 (corresponding to 
Sr ~ 0.04 ... 3, where r is the turnover rate of the ran- 
domly forced velocity fluctuations), we have found that the 
dynamo growth rate is 7 ex S* and the characteristic length 
scale of the generated mean magnetic field is Ib cx I/a/S'. 

The first key result of this pap er, compared with the ear- 
lier study by lYousef et al.l (l2008b . is that the shear dynamo 
works both in the nonrotating case and for the case of Ke- 
plerian rotation (Q = —2S/3). There does not appear to be 
much difference, qualitative or quantitative, between the ro- 
tating and nonrotating cases, although perhaps it would be 
interesting to look at non-Keplerian cases and try to identify 
the role of rotation via a parameter scan in Q independent 
of the one in S. 

The second key result, claimed on the basis of only 
a preliminary study, is that the shear dynamo works both 
for situations that are sub- and supercritical with respect to 
the fluctuation (small-scale) dynamo. In the latter case, the 
overall magnetic energy grows very quickly due to the fluc- 
tuation dynamo effect independent of the presence of the 
shear. Imposing the shear on the magnetohydrodynamic tur- 
bulence resulting from the saturation of the fluctuation dy- 
namo, leads to the emergence of magnetic fields that have 
spatial scales larger than the outer scale of the turbulence 
and that fluctuate on very long time scales compared to the 
turbulent turnover time. 

Thus, the shear dynamo effect appears to be quite gen- 
eral and robust. As the combination of a shear flow and tur- 
bulence is a very common feature in astrophysical systems. 
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the shear dynamo potentially represents a generic mecha- 
nism for making large-scale fields. While much needs to 
be understood about its properties before its relevance to 
real astrophysical systems can be more than an appealing 
speculation, the simplicity of the idea of the shear dynamo 
certainly makes it a worthwhile object of study. It is also 
important to determine how generic the shear dynamo is 
and how it combines with other large-scale features present 
in real astrophysical systems: various differential rotation 
laws, temperature and density gradients, linear instabilities, 
etc.. 

Studies in this vein are already being under- 
taken. For example, recent numerical experiments by 
Kapyla. Korpi & BrandenburgI (12008^ have shown that 
large-scale dynamo action is also possible in local simula- 
tions of magnetoconvection with imposed horizontal shear. 
They also report a growth rate j ^ S and find large-scale 
dynamo action for magnetic Reynolds numbers above the 
critical threshold for t he fluctuation dynam o. Another topi- 
cal recent study is by lOressel et al.l (12008^ . who simulated 
the supernova-driven galactic turbulence and found that 
they needed to impose a linear velocity shear to obtain the 
amplification of a large-scale field. These studies clearly 
demonstrate the key role of shear in producing a mean-field 
dynamo. However, in comparing their results to ours, one 
has to keep in mind that their simulations had rotation and 
vertical stratification, so the turbulence in these simulations 
is likely to be helical and may also host an a effect. 

In motivating our choice of Keplerian rotation law, 
we mentioned the possible relevance to accretion-disk tur- 
bulence, which is believed to b e driven by t he mag- 



netorotational instability (MRI) (IBalbusI l2003l) . Given 
a (weak) large-scale field, the MRI will generate ve- 
locity and magnetic-field fluctuations at small scales. 
These fluctuations, in conjunction with Keplerian rota- 
tion and shear, must then amplify the large-scale field 
to close the loop. The mechanisms for such an "MRI 
dynamo" have been discussed and simulated for some 
time (e.g.. iBrandenburg et al.l 1 19951: iFromang et al.l 2007 : 



Hawlev. Gammie & Balbusni996|: iLesur & Qgilvij 12008: 



Rincon. Ogilvie & Proctoj 120071: bincon et al.ll2008l) . It is 
tempting to observe in the context of the results reported 
above that a combination of small-scale turbulence (mag- 
netohydrodynamic turbulence in the case of the MRI) and 
large-scale shear does indeed appear to work as a dy- 
namo giving rise to a large-scale azimuthal magnetic field 
(By in the shea ring sheet model). We note, however, that 
Lesur & Ogilvie (2008), who have analysed this process 
in detail, find some important differences between what 
happens in MRI-driven shearing sheet simulations and the 
forced case studied by us. 

To conclude, we believe that the discovery of shear dy- 
namo has opened a number of new and exciting avenues of 
research and produced some promising leads towards un- 
ravelling the ways in which cosmic magnetic fields emerge. 
Further investigations will help assess the range of appli- 



cability and relevance of the shear dynamo effect and the 
physical mechanisms that are responsible for it. 
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